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N-Time-Gate Data-Type for TPSF-Based Optical Imaging 

TECHNICAL FIELD 

This invention relates to the field of temporal point 
spread function (TPSF) based imaging in which objects which 
diffuse light, such as some human body tissues, are imaged 
using signals resulting from the injection of light into 
the object and detection of the diffusion of the light in 
the object at a number of positions while gathering TPSF- 
based data to obtain information beyond simple attenuation 
such 3S scattering and absorption. More particularly, the 
present invention relates to a specific data-type, N-tlme- 
gate, which can be obtained from TFSF-based data. 

BACKGROUND OF THE INVENTION 

15 Time-domain optical medical imaging shows great promise as 
a technique for imaging breast tissue, as well as the brain 
and other body carts. The principle of the technique lies 
xn the analysis of the temporal point, spread function 
(TPSF) of an injected pulse of light which is diffused in 

20 an object of interest (001), and in • the extraction of 
informaticn from the .TPSF to construct a medicaid useful 
image.. The (complete) TPSF is defined as the light ■ 
.intensity of a given wavelength presented as a function of 
time t (t e st«>») that can be observed at a given spatial 

25 location after a very short pulse of light produced at 
another spatial location at the time t 0 . 

Specific types of data extracted from the TPSF are known as 
'Data-types'. For example, 'early time-gated attenuation' 
is a Data-type that can be extracted from the TPSF which 
improves the image spatial resolution over continuous wave 
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methods. However, it is unclear whether such improvement, 
in mage spatial resolution are adequate for diagnosing 
certain diseases, such as breast . cancer , based \ Qn 
morphology . 

An alternative approach is to use the TPSF to decouple the 
light attenuation into absorption and scattering 
components. This extra information, which generally cannot 
be obtained with continuous wave methods, may be clinic y 
useful, in order to achieve this it is customary to extract 
appropriate Data-types from the TPSF. Researchers, in the 
time-domain field, have looked at data-types such as the 
meantime and higher moments of the TPSF (Arridge, 1939) 
Whilst .data-types for optical data acquired in the 
frequency domain have included phase-shift and amplitude. 

There are several techniques for measuring the TPSF arista 
from a light pulse that has traveled through an OOI such as 
breast tissue. Certainly the most information one can 
acquire is by measurxng a complete TPSF. This information 
can also be acquired in the frequency domain, although 
current hardware limitations mean that time-domain hardware 
is capable of acquiring information over a larger bandwidth 
thar: its frequency domain counterpart. 

However, acquiring a practical approximation /(,), / e[a,b] to 

i 

the complete TPSF, such that //<**.]/<,)*, has implicaticns 

in terms of long acquisition times for clinical systems, 
ana possibly acceptably long data-processing for image 
reconstruction algorithms. 

I. would therefore be desirable to provide methods 
overcoming the limitations of the prior art. 
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SUMMARY OF THE INVENTION 

For the purpose of image construction, the entire TPSF 
derived from a pulse of light can be used.- The streak- 
camera has commonly been used for this purpose. However- the 
detection, surface of a streak camera may be limited. Time- 
gated cameras such as intensified charge coupled device 
(ICCDj typically have greater .detection surface areas but 
acquisition of complete TPSFs using such camera may be time 
consuming. 



It is therefore an object of . che present invention to 
provide a method and apparatus for TPSF-based imagine of an 
object in which acquisition times are reduced while still 
extracting a sufficient amount of desired information from 
the TPSF of a collected optical signal to be used in 
15 producing an optical image. 

According to a broad aspect of the invention this objective 
is achieved by detecting selected intervals, or N-time- 
gates, of a TPSF. The detection of the selected time-gates 
xs effected, using time-gated detection techniques and 
20 cameras. 

In a further aspect of the invention, che acquisition time 
is. substantially reduced by simultaneously detecting 
selected time gates of a TPSF. 
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Thus in one aspect of the present invention there is 
provided a method for simultaneously detecting N-time-gates 
Data-type of a TPSF of a collected optical sicnal, the 
method comprising pulsing the object with a light pulse 
delivered at one or more injection ports, coll.ctino a 
plurality of optical signals based TPSF'. at collection 
ports, introducing delays in propagation of r.he optical 
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signals to produce staggered TPS~. S ari H 

yy ea ie5i s an ^ simultaneously 
detecting selected time-gates. 

•in an embodiment of the invention the delavs in the 
propagation of the optical signals are pmi ^ fay ^ 
b optical fibers' of different lengthg tQ * 
Signals. 

BRIEF DESC RIPTION OF THE DRAWINGS 

Further features and advantages of the present invention 

10 r°" !e aPParent f — ' the follows detailed 

10 description, taken in combination u1 » h -v. 

°- nat10 " Wlth the apoended 

drawings, in which: 

Fig- 1 is a diagram showing an embodiment of tr= s yst en- 
used to generate a pulse of li ght 3nd t0 produce a ^ 
after injection of the light in an object; 

Fig- 2 is a graphical representation, of a TPS F showing 
arbitrarily defined time-gates; 
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Fig. 3 is an example of a system that h. 

^iiau c<*n oe usee *■ o 

construct an optical ina 9 e according to methods of the 
present invention; 

20 4 13 3 8CheM tiC -P-sentation of an embodiment of 

--he system used to simultaneously detect two or more tire- 
gates; 

XPsrs! " * ?raPhiCai r ' Pre " n "-- i - of three starred 
.-•«. a IP^F. .. on , with a first derivative; and 



Fig. 6B is a plot of the first derivative of Fig. 6A along 
with a second derivative. 

It .will be noted that throughout the appended drawings, 
like features are identified by like reference, numerals. 

DETAILE D DESCRIPTION OF THE PREFERRED EMBODIMENT 

In accordance with the present invention, there is provided 
•a method and apparatus using a w-time-gate Data-Type for 
TPSF-based optical imaging. 

Figure 1 schematically represents the events giving rise to 
a T?SF. A pulse of light 2 is generated, at a source 10 and 
injected in an OOI 4 at injection port 14. The OOI may be a 
body part such as a breast but it will be appreciated that 
the method and apparatus described herein can be applied 1 to 
any suitable object. The light interacts with the medium by 
being absorbed and/or . scattered and exits the OOI with' an 
attenuated intensity.. The exiting attenuated light can be 
collected at collection port 18 as a function of time to- 
produce the TPSF 6. The shape of. the T?SF depends on the 
type of interactions that occurred between the light and 
the medium. Images (maps of optical properties) may be 
•constructed by extracting appropriate Data-types from the 
TPSF and by inputting them to an appropriate image 
reconstruction algorithm, i.e. one which forward models the 
daca and solves the appropriate inverse problem. 

In one embodiment, the present invention provides a method 
for detecting only selected intervals of a TPSF. In Figure 
2 it can be seen that a TPSF can be subdivided in n 
arbitrary intervals 8 which will be referred to as time- 
gates (eg-; hereinafter. The time-gate is defined as an 
interval [a,i>)<z (r 0 ,eo). where a and b are predefined constants. 
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Therefore one can consider a complete TPSr as being 
constructed from many individual time-gate signals. The 
actual time points that define the beginning and the end of 
a time-gate are preferably determined relative to the time,. 
5 t 0f at which the light pulse is injected into, the object. By 
detecting the signal of sequential time-gates using time 
ga^ed detection techniques and cameras such as, but net 
limited to, a time-gated intensified charge coupled device 
.• (ICCD), it . is possible, to acquire the complete TPSF. 
10 However, acquisition of a complete TPSF in this manner 
requires long acquisition times which will be proportional 
to the number of time-gates acquired. In one aspect of the 
present invention the image may be constructed by acquiring 
only certain time-gates. As a result the acquisition time 
is greatly reduced. Overlapping time-gates can be allowed 
if the forward model is capable of taking this into- 



account 



The selection of the time-gates relative to t c that will be' 
used to. construct ar. image may be based on the 
characteristics of the object such as the thickness and the 
anticipated optical properties. In the absence of prior 
knowledge about the anticipated optical properties of the 
object, a practical approximation to the complete TPSF 
could be obtained and be used for future selection of the 
25 time -gates . 

Ir. a preferred embodiment certain specific time gates or 
time points (infinitely small time gates) of the TPSF at 
which TO data will be acquired can be selected based on the 
following approach: 



30 Based on the physics of diffuse light propagation, 1; 

be assumed that a noiseless TPSF is a continuous and 



can 



infinitely differential, function defined on r\ Critical 
points of a function can be identified via zeros of .his 
function's derivative. Due to the fact that a TPS? has a 
characteristic shape (see 6 Figure 1), the ,-th derivative 
will yield a set of points where it became zero, which 
corresponds to extreme values of the underlying- function. 
The first subset ,,={,-,] corresponding to the 1"' derivative 
describes the unique TPSF xaximum. The second subset 
^-K'i f wl corresponding to the 2^ derivative reveals two 
points where the curvature directionality changed 
Derivatives of the higher order yield information on riore 
subtle TPSF shape features. The subset ,] corresponding to 
the i-th derivative can be defined as s, = {*,,;.. ;,,..} . 

Given a TPSF /(,) represented with a set of M values ' <K»n) , 
f-ifOlflfA. -.At*)), ,,*«,- t M <b, one can perform 
recursive numerical "differentiation" of the TPSF and of 
its' resulting "derivatives" by computing the quotient: 

0 . /(',»,)-/(/..,) 

which would approximate the value cf the derivative /'(/.). 

When plotted, characteristic sequence of functions will 
emerge as shown in Figures 6A and 68. One observes that in 
this example s,-{3:.5), s*={20.5; 4C.5), Sj ={14.5; 25.5; 
48.5), the values being the time point indices rather than 
absolute time. 

in order to obtain an adequate TPSF representation with n 
points (or time gates), one would want to reflect the TPSF 
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shape as fully possible. The seguence {,,;,,,..} ilWs £or 
■ id.ntxfic.Won of the «o, t ' „ lmM points in t „. ^ 

drmntshing significance es the subset nuafcer increases ' 
• KOCe elKMat » UUhin « ' *v«, subset « have e qua .' . 

5 Significance with e noiseless T ?SF . However, elements of , 
■•y be ordered by assigning higher priority to elements 

that have higher maxk _c I <~ ^ 

y ^1 ln °rder to minimize the 

correlation cf experimental errors. 

One can consider the- midpoints cf a set of time-gates 
i.e. ^ = k-0.5r />I ,, M+ o.5,J..4 J/i -o.5,,,,,, +0 .5,,.)} / where r<j is the . 
width of a respective time gate. 

The procedure outlined above can yield S' -{r,';^/-! K to- 

a given IPSF having index i. m order to reflect the most 
relevant features of a set of A* TPSFs, one must arrive to 
S^s^^s? first, to ^^^...^ second and SQ cn 
The set {S,;S 3 ,...} provides a sequence whose « sta-ting 
elements should be used to represent the set of -TPSFs in -an 
optimal manner, 

As mentioned above, s. might be replaced with a s-t o- 

20 intervals m order to analyse time-gates instead of w„e 
points. 

in a preferred embodiment of the invention the steps to 
select the appropriate time gates are as follows: ' 

- Acquire practical aporox^at-on? 

ayr'io^. 1 ..a:.ons ..o t;-»e complete TPSf* 

^ zor all source-detector pairs that are capable of 
functioning simultaneously. At least one typical source 
position through the bulk medium , heuid be used. Th . ' 
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robustness of the subsequent analysis increases, as 
typica. source positions are included. 

■ - Approximate all measured TPSFs either based on a 
relatively simple model of photon propagation, e . g 

5 asking that propagation took place in a homoaeneous 
medxum, or interpolate the measured TPSFs to ensure the ' 
continuity and differentiability of the result. 

- Perform analysis yielding the set as described 

above in order to determine the optimal n time-gates. 
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n is an 
via 
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20 
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- Acquire data with the determined n time-ga^es 
input parameter, m practice it can be determined 
estimation of the resulting data set size rhat can be 
utilized to reconstruct images within a desired time frame. 

Temporal resolution r 0 , i.e. the minimal distance between ,,. 
and t,„ that yields virtually "uncorreiated" values /(,,) and 
/(0, is determined by the limitations of the data 
acquisition Hardware and is, therefore, an input parameter 
It x. preferable to perform initial data acquisition with 
time-gate width and time steps of not more than 0.5*r e 
Eventual width ri >r 0 of the time-gate having index i is 
automatically chosen to . ensure capturing the most relevant 
TPSF features revealed by the analysis. 

The above described method for selecting time oates may be 
usee i„ the context of a broader process for select 
appropriate time-gates when xmaging a particular object Ir 
one embodiment this broader process may comer i M -r- use o- 
banks of TPSF-s obtained from a plurality of objects having 
given recorded characteristics and for which the derivative 
method described above has. been applied to select 
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appropriate tine-gates. Thus, the selection of appropriate 
time-gates for a new object being imaged may then be made 
by matching certain predetermined characteristics of r h ~ 
object with those recorded in the data banks. When a "good" 
5 match has been found the time-gates that have been 
identified for the corresponding object' in the data ban, 
are used to acquire the TPSF of the object being imaga d. 
For example, in the case of breast imaging, TPSF's can be 
acquired from a plurality of patients having different 
10 characteristics that can be recorded such as age, breas^ 
density, disease state etc. The derivative method is then ' 
applied to. these TFSF's to select appropriate time-gates 
and the identity of the time-gates are associated with each 
patient. When it is desired to image the breast of a new 
15 patient the appropriate time-gates may be selected by 
simply matching the characteristics of this patient with 
those m the data bank thereby effecting the selection of 
the time-gates. it will be appreciated that the 
characteristics used in the data bank maybe chosen to 
20 reflect their relevancy to optical imaging. 

An embodiment of the system that can be used to perform the 
data acquisition using the method of the present invention 
is. schematically represented in Figure 3. The oulsed light 
source 10 has an output optically coupled via a swirch 12 
to one or a plurality of injection ports 14 of a support. 
The injections port are preferably positioned at a number 
of fixed positions over the imaging area aithouah the 
inaction ports may alternatively be moveable over the body 
surface, provided that the body part 16 beina imaoed is 
30 immobilized. As is known in the art the injectio'n and 
collection ports may directly contact the body or a 
coupling medium may be used between the- body and th- 
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section/collection pets, m Figure 3, the collection 
porta 18 and support are arranged in transmission mode -for 
breast imaging. However, it is al so possible to have " an 
in 3 ection ports/obj ect/collection ports geometry such that 
the., collection ports are on the same surface as the 
injection ports. The light source may be a polychromatic 
source or a monochromatic source such as a laser The 
pulses preferably have a duration of about 1 to 100 
picoseconds and an average power of loo mW. ' 

The collected optical signals may be communicated by one or 
mere waveguides 20 such as 400/440 micron graded index 

mul timod«s optical f i*Aav. c *- ^. ~~ _ 

optical .iseri to one or more detection positions 

of a time gated detector 22 such as a gated ICCD camera 

ror example a PicoStar" camera by LaVision. 

A large number of pulses nay be. injected and their 
corresponding camera signals detected and processed bv 
imaging computer 24 to determine the value of at least one 
time gate and, in the case of simultaneous detection, two 
or more selected time gates of a TFSF. Such data are 
gathered for . sufficiently large number of detecting oor, 
positions to ensure robust reconstruction of an imag- o* 
tne object. The resulting image can be disolayed on a 
monitor 26 or a hard copy can be obtained at 28 . 

It will be appreciated that the relative positions of the 
injectors and detectors may be modified as needed to 
optimize the quality of the image. Thus different 
geometries cf the injection ports/collection ports/object 
assembly may oe used, for example by placing, collection 
ports either proximally or distally or a combination of 
botn, depending on several factors such as, the thickness 
of the object, the nature of the medium comprising ^be 
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generated from the data collected according to the p , Mlrt 
invention win benefit from a plurality of source detector 
locations with respect to the position of the pixel or 
voxel, as well as a plurality of time-gates. 

in an aspect of the present invention the acquisition r llPa 
can be further reduced by simultaneously detecting a 
assired number of time-gates of one or more TPSF's Thi* 
can be achieved by collecting the light signal exiting from 
the OOI at several {m) locations. Temporal delays are t-h-r 
introduced in the propagation of these optical signals such 
that the ■» TPSc's reach the time-gated detector ir a 
staggered manner. Therefore, using a time-gated camera, it 
is possible to simultaneously detect m different tire-ga^e. 
15 cf the TPSF's. ' 

This embodiment of the invention will be better understood 
with reference to Figures 4 and 5. In . Figure 4 a schematic 
diagram of a possible arrangement of a system used to 
simultaneously detect » time gates is deoictec. In this 
embodiment optical fibers 30 . are used to ooilect and 
channel the light exiting the 001 4 to a time-gated 
detector 32. The fibers are of different lengths and 
therefore the time required for the optical sional to 
travel from the collecting point to the detection point at 
25 the detector is different for each fiber thus introducing 
decays m the propagation of the optical signal. Ea-T 
fiber, or bundle of fibers, is coupled to a dlffer#nt 
detection position 34 at the detector. Consequently, the 
tune required for a particular time gate, tg , to reach th, 
detection position ■ in the detector will depend on the 
length of the fiber through which it is propagated. Figure 
5 ~s a graphic representation cf three staggered TPSF's 
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propagated by three different optical fibers n, F2 , and F3 
• having different lengths with length of F1 > length of P2> 
.ength of „. It can be appreciated that a given time-gate, 
say tc3, reaches the detector at a- different time deoendl^ 
on wnether ic is propagated by F1 , . F2 or F3 ; Now \ f ^ 
detector is time- g ated to detect only the signal during the 
interval dexxned by tl and t2, different time-gates wUl 
oe defected f or each fiber, in our example t{7 3, tg9 , a " nd 
tg7 are sxmuxtaneously detected. x, . can be ssen th „ 
number of time gates that can be detected simultaneously is 
equal to the number of fibers. The fibers can be bunded 
wxthxn a defined area of a given detector and ^ 
desired time-gates simultaneously with one 'snapshot Q l 
the tim.-g.ted camera. This leads to a decrease in the 
acquisition time by a factor of approximately n compared to 
serially acquired time-gates. 

It will be appreciated that the length of the fibers can be 
selected to detect the desired time gates (or the desired 
sections of the TPSF) . To refer to our example, if one 
wanted to simultaneously detect tg7, tg5 and , the 
length of F1 COald be increa£ed accordingly . ^ 

will axso be appreciated that the simultaneous detection of 
the desired time-gates can be effected by directing each 
nber to a different detector and by adjusting the time 
gate of each detector so as to detect one of the aesired 
time gate. Typically, an increase of 11 cm in f iber • lenat , 
will induce a time delay of about 500 ps. 

T dei3yS in P-P^ation of the optical signals may 

also be introduced by using variable delay" optical 
waveguides that are well know in the art. 
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It wili be appreciated thai: the signal to noise ratio (SNR) 
may lir.it -he number of fibers and thus -he number of time 
gates that can be simultaneously detected, in practice, 4t 
is a decision between how many detector positions are 
.required, the number ■ of time-gates simultaneously required 
and the total number of fibers which can be packed onto the. 
gated camera. Referring back to Figure ■ 4, it will be 
further appreciated that if the fibers or the fiber bundles 
35 are sufficiently close- together, the TPSF's may be 
substantially identical. 

The embodiment (s) of the invention described above is (are) 
intended to be exemplary only. The scope of the invention 
is therefore intended to be limited solely by the scope- of 
the appended clairr.3. 



